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Abstract 

Variation of constants in the very early universe can generate infla- 
tion. We consider a scenario where the strong coupling constant was 
changing in time and where the gluon condensate underwent a phase 
transition ending the inflation. 
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1 Introduction 

Alternatives to inflationary cosmology [U E] include varying speed of light 
(VSL) theories [3J- Usually all inflationary models are based on using new 
fundamental scalar fields, the 'irmatons', whose nature is still unknown. Some 
models change the matter content of the universe, while others give the in- 
flaton geometrical interpretations within brane settings [3J. VSL scenarios 
may solve the cosmological problems usually tackled by inflation ( "horizon" , 
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"flatness" and "structure formation" problems) without introducing infla- 
tons, whereas many inflationary models lead to variation of other constants 
of nature [5]. In this paper we will follow the reverse route that variation 
of constants in the very early universe can generate inflation. In an earlier 
work [6 J , we considered a Bekenstein-like model for the QCD strong coupling 
constant as introducing a scalar field e expressing the time variation of as- 
We found that experimental constraints going backward till quasar forma- 
tion times rule out as variability. However, when this model is implemented 
in the very early universe, the scalar field e can play the role of an infla- 
ton, and one can realize a consistent inflation scenario with suitable value 
for the gluon condensate. We find that the 'time varying' QCD lagrangian 
leads naturally to a monomial quadratic potential like the chaotic scenario. 
However, while the large values of the inflaton matter field plague the latter 
scenario, they just amount in our model to a reduction of the strong charge 
by around 10 times during the inflation. An exit way can be achieved if the 
gluon condensate suffers a phase transition reducing its value to its current 
value ending thus the inflation. We will not dwell on the possible mecha- 
nisms for such a phase transition to occur, but wish to concentrate on the 
conditions our model should satisfy in order to present a consistent set up 
able to accommodate the recent measurements from the Cosmic Microwave 
Background (CMB) E] and the WMAP results PE)]. 

2 Analysis 

We follow the notations of [6]. Our starting point is the 'time varying' QCD 
Lagrangian 

Lqcd = L e + L g + L m 

= -^^-\Tr{G^G, v )+^ U \^D,-M f )^ (1) 

where I is the Bekenstein scale length, e(x) is a scalar gauge-invariant and 
dimensionless field, with the 'variable' QCD coupling given by g(x) = goe(x). 
The covariant derivative is = 8^ — igoe(x)A^ and the gluon tensor field is 
G% = \[d,{eAl) - d v {eAD + g^r hc A\Al\ 

Assuming homogeneity and isotropy for an expanding universe we con- 
sider only temporal variations for as = = as e 2 (t). One gets the follow- 
ing equations of motion 

(^);M-^o/ abc GT^ + E^ a 7^ = (2) 
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where a(t) is the expansion scale factor in the R-W metric. 

Subtracting the total derivative AT a/3 = d v (— ^-eA al3 ^J from the canon- 
ical energy-momentum tensor Q^pjd 13 ^ — g a ^L we get the gauge-invariant 
energy momentum tensor 
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In the analysis below, we consider only the gauge fields and the scalar 
field, i.e. we drop the matter fields in the radiation dominated very early 
era. We decompose our energy momentum tensor into two parts: the gauge 
part and the e-scalar field part 
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Here all the operators are supposed to be renormalized and it is essential 
in the inflationary paradigm that quantum effects are small in order to get 
small fluctuations in the CMB. 

The contribution of the scalar field to the energy density p e = T e and to 
the pressure = g^p e are 



Pe 



(6) 



On the other hand, the gauge field contribution can be decomposed 
into traceless and trace parts. 



Pg = Pl + Pg ( 7 ) 
Pg = fg+P^ (8) 

where p r g ,p r g are the density and the pressure corresponding to the "trace- 
less" part of the gauge field satisfying p r g = 3p r g , while the trace part of the 
gauge field energy-momentum tensor is proportional to g a p and behaves like 
a 'cosmo logical constant' term: 

P T 9 = ~P T 9 (9) 
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This equation is reminiscent of 'ordinary' inflationary models. However, to 
compute the trace part of the density one needs a 'trace anomaly' relation for 
our 'time varying' QCD. Since the energy-momentum tensor Tfg is identical 
in form to "ordinary" QCD and since the trace anomaly which involves only 
gauge invariant quantities should, by dimensional analysis, be proportional 
to G 2 , we take it to be the same as in "ordinary" QCD (we have checked 
that changing the numerical value of proportionality will not alter the con- 
clusions). Thus we take, up to leading order in the time- varying coupling 
"constant" as = as e 2 , the relation [TT] : 

?7 = P 9 ~ 3P 9 = G% (10) 

This leads to: 

p r = _|^ f2 < G , > (n) 

As we said before, equation suggests, in analogy to ordinary inflationary 
models, that the QCD trace anomaly could generate the inflation. For this, 
let us assume that the "trace-anomaly" energy mass density contribution is 
much larger than the other densities: 

p T g»Pe,p r g P~P T 9 (12) 

Then, equation ( TTTT) tells that the vacuum gluon condensate < G 2 > should 
have a negative value which is not unreasonable since the inflationary vacuum 
has "strange" properties. In ordinary inflationary models, it is filled with 
repulsive-gravity matter turning gravity on its head [TJ]. This reversal of 
the vacuum properties is reflected, in our model, by a reversal of sign for the 
vacuum gluon condensate. 

Now we seek a consistent inflationary solution to the FRW equations in 
a flat space-time: 



, a.. 2 8ttG]\[ 
V ~ 3 



P (13) 
(p + 3p) (14) 
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where Gn is Newton's constant. The first FRW equation with ( TTT]) will give 



H =- = eJ^G N \<G 2 >\ (15) 
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On the other hand, the equation of motion (eqJH]) of the scalar field can be 
expressed in the following way: 



i +aH i-fy = *±£> (16) 

This equation differs from the ordinary 'matter' inflationary scenarios in the 
term (^) 2 . However, for "slow roll" solutions we neglect the terms involving 
- e and (^) 2 to get 

3He = l ^f^e =-V\e) (17) 

which is the same as the "slow roll" equation of motion of the inflaton in 
ordinary scenarios. In our model, the "slow roll" condition can be written 
as: 

5=\4r\ = qMt-)H «1 (18) 

We set e/, the value of e-field at the end of inflation tf, to 1 so that the 
time evolution of the strong coupling terminates with the end of inflation 
and we expect for "slow roll" solutions that Cj, the value of e at the initial 
time of inflation ti corresponding to when the CMB modes freezed out, to 
be of order 1. If the gluon condensate value < G 2 > stays approximately 
constant during much of the inflation, the changes of the Hubble constant 
and the energy mass density are not very large. In this case we have 



*) = ^ 337^(i) 2 <l< G2 >l 1/2 ('-'<) < 19 > 
and , as in chaotic scenarios, we get a simple quadratic potential: 

V{e) = >l e2 (20 ) 

One can make explicit the correspondence between our model with e- 
scalar field and the chaotic scenario with an 0-inflaton matter field. Com- 
paring equations ( |T5l) and (1T7I) with the corresponding relations in ordinary 
inflationary models: 

3H<j) = -V'(0) (22) 
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we find the relations between (e, V(e)) and (4>, V((f>)): 



with 



V(0) 
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(23) 
(24) 



3 Results and Conclusion 

Now, we check that our model is able to fix the usual problems of the standard 
(big bang) cosmology. First, in order to solve the "horizon" and "flatness" 
problems we need an inflation ^jA- of order 10 28 implying an inflation period 
At = t f — t i such that 

H At ~ 65 (25) 

Furthermore, it should satisfy the constraint 

10" 40 s < At < l(T 10 s (26) 

so that not to conflict with the explanation of the baryon number and not to 
create too large density fluctuations [131 E]. The bound 10 _10 s corresponds 
to the time, after the big bang, when the electroweak symmetry breaking 
took place. Presumably, our inflation should have ended far before this time. 
Thus, from equations ( |25l) . ( 1261) and ( fl5j) we get the following bounds on 
| < G 2 > |: 

3 x 10 7 GeV 2 < e| < G 2 > \ 1/2 < 3 x 10 37 GeV 2 (27) 

In order to determine e^, we have = M- ~ — ~ m'*v' w hich gives, 

using equations (fill and (fTTl) . the relation: 

65 ~ln^f= (^) 2 ^; 2 -D (28) 

Next, comes the "formation of structure" problem and we require the frac- 
tional density fluctuations at the end of inflation to be of the order 
10 -5 so that quantum fluctuations in the de Sitter phase of the inflationary 
universe form the source of perturbations providing the seeds for galaxy for- 
mation and in order to agree with the CMB anisotropy limits. Within the 
relativistic theory of cosmological perturbations [T5], the above fractional 
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density fluctuations represent (to linear order) a gauge-invariant quantity 
and satisfy the equation 

— I — I 1 

p 

where 5M represent the mass perturbations and where the initial fluctuations 
are generated quantum mechanically and are given by: (T5J 03] 

5M _ V(*)H _ V'(e)H 

w k ~^p~ - ^~ip~ m 

whence 

* , 6M , _ , V'H ,1 

10 " l ¥l" = ^1 — I^^Trtk (31) 

In order to evaluate (p we use the energy conservation equation: 

P + 3(p + p)- = (32) 
a 

and after substituting p ~ p J we get 

In fact, the energy conservation equation can be used to solve for p r g and we 
could check that 

fr g (fr g )~pM ~ Sxp^SxpT) (34) 

where 5 = |-^| ~ ^(r^) 2 an d so ) when the "slow roll" condition fj 1 8p is 
satisfied, our solution assuming the predominance of the "trace-anomaly" 
energy mass density is self-consistent. Substituting equation fl3"3l in fl3T|) and 
using equation ( TTT1) we get 

^ ~ 9^a 5o 10V|<G 2 >|^^ (35) 

Hence, taking Gn ~ lCr 38 Gey~ 2 we obtain 

I \<G 2 > 



2 2 



L P W^GeV 2 



(36) 
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and combining this last result with (1271) . we get 



1(T 27 « --j- < 10 3 (37) 
e L„ 



The "slow roll" condition (1181) is consistent with the upper bound, while the 
lower bound restricts q not to be too large. 

On the other hand, it is possible to calculate the spectral index of the 
primordial power spectrum for a quadratic potential as follows: 



Ml,V" 2 . I N2 1 
n - 1 = -4?y where r? = = h— — = d 3£ 



and we find: 



n = l-—{-^-) 2 \ (39) 
ixy L p e\ 

The inflation would end (ey = 1) when the "slow roll" parameter r] — 8 — 

1. We should evaluate the QCD coupling constant ots (fA = ^-2 — at 

AMtI — ) 

QCD 

an energy scale corresponding to the inflationary period. We take this to 
be around the GUT scale ~ 10 15 GeV and (3q = 11 — |n/ = 7 (the weak 
logarithmic dependence would assure the same order of magnitude for as 
calculated at other larger scales). With Aq C d ~ 0.2GeV [T7] we estimate 



a s ~ 0.025, and so we get 



;-^) 2 ~ 10- 1 (40) 

LJ 



P 



This is in disagreement with Bekenstein assumption that Lp is the shortest 
length scale in any physical theory. However, it should be noted that Beck- 
enstein's framework is very similar to the dilatonic sector of string theory 
and it has been pointed out in the context of string theories [T8] that there 
is no need for a universal relation between the Planck and the string scale. 
Furthermore, determining the order of magnitude of jr- is interesting in the 
context of these theories. 

From (1281) . we have e; ~ 11, and then using ( |39|) . we have n = 0.97 which 
is within the range of WMAP results [HI [10] . 

The model reproduces the results of the chaotic inflationary scenario. 
However, the shape of the potential was not put by hand, rather a gauge 
theory with a changing coupling constant led naturally to it. Moreover, 
in typical chaotic models, the inflaton field starts from very large values 
(<pi ~ 15Mpi) and ends at around IMpi. One might suspect whether field 
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theory is reliable at such high energies. Nonetheless, this problem is absent 
in our model since the large values have another meaning in that they just 
refer to a reduction of the strong coupling by around 10 times during the 
inflation. 

Furthermore, chaotic inflations get a typical reheating of order T r h ~ 
10 15 GeV, and one might need to worry about the relic problem. Similarly, 
equation (I3BT) leads in our model to a gluon condensate | < G 2 > |j ~ 
10 62 GeV 4 at the start of inflation. From equation (TT9l) . we see that this 
corresponds to an inflation time interval At ~ 10 _35 s satisfying the con- 
straint ( 1261) . If the gluon condensate stays constant, as we assumed in our 
analysis, we will have the same reheating temperature as in chaotic models 
(T rh ~ p(ty) 1//4 ). However, we should compare this value for < G 2 > with its 
present value renormalized at GUT scale ~ 10 15 GeV which can be calculated 
knowing its value at lGeV [6] and that the anomalous dimension of a>sG 2 is 
identical ly zero. We get 

< G 2 (now, 11 ~ 10 15 GeV) >~ lGeV A (41) 

which represents a decrease of 62 orders of magnitude. 

This can give a possible picture for an exit scenario. Lacking a clear theory 
for the non-perturbative dynamics of the gluon condensate, we consider its 
value | < G 2 > | depending on energy, and thus implicitely on cosmological 
time, as given by the standard RGE which turns it off logarithmically at 
high energy. However, we can furthermore assume the condensate value to 
depend explicitely on time during inflation: 

<G 2 {E,t) >=< G 2 (E(t)) > f(t) 

where < Gq(E) > is the piece determined by the RGE, the unknown function 
f(t) should be such that it varies slowly during most of the inflationary era, to 
conform with an approximately constant huge and negative value of < G 2 >, 
while at the end of inflation it causes a drastic drop of the condensate value 
< G 2 > to around zero. The energy release of this helps in reheating the 
universe, while reaching the value leads to a minute "trace-anomaly" energy 
mass density (equation [TT]) ending, thus, the inflation. The other types 
of energy density would contribute to give the gluon condensate its 'small' 
positive value of (HTi) . and the subsequent evolution is just the standard one 
given by RGE. Surely, this phenomenological description needs to be tested 
and expanded into a theory where the concept of symmetry breaking of such 
a phase transition for the condensate < G 2 > provides the physical basis for 
ending the inflation. Nonetheless, with a test function of the form f(t) = 
— /5 2 tanh 2 (e — 1) with f3 ~ 10 31 , one can integrate analytically the equation 



9 



of motion, and in "slow roll" regime we have e = 1 + Arcsinh(exp[— a/3t]) 
with a(3 ~ 10 12 GeV. The graph in Fig. [T] shows the time evolutions of 
the condensate and the e-field, which agree with the required features. This 
example is meant to be just a proof of existence of such functions, and the 
temporal dependence of the condensate could be of complete different shape 
while the whole picture is still self-consistent. The issue demands a detailed 
study for the condensate within an underlying theory and we do not further 
it here. We hope this work will stimulate interest in the subject. 



Figure 1: Temporal evolution of the condensate | < G 2 > | (thick line) and 
the e-field (thin line), for the choice f(t) = -/3 2 tanh 2 (e - 1). The < G 2 > 
scale has been adapted so that to visualize both graphs together. 



Acknowledgements 

This work was supported in part by CONICET, Argentina. N. C. recognizes 
economic support from TWAS. 

References 

[1] A. H. Guth, Phys. Rev. D 23, 2, 347 (1981) 

A. Albrecht and P. Steindhardt, Phys. Rev. Lett. 48, 1220 (1982) 

[2] A. Linde, Phys. Lett. B 129, 177 (1983) 

A. Linde, D. Linde and A. Mezhlumian, Phys. Rev. D 49, 4, 1783 
(1994) 




time 



t_i 



t_f 



10 



[3] J. Maguejo, Rep. Prog. Phys.66 (2003) 2025 

J. Moffat, Int. J. P. D 2, 3, 351 (1993); Foundations of Physics, Vol. 
23, 411 (1993) 

A. Albrecht and J. Magueijo, Phys. Rev. D 59, 04, 3516 (1999) 

[4] J. Khoury, B. Ovrut, N. Seiberg, P. Steindhardt and N. Turok, Phys. 
Rev. D 65, 086007 (2002) 

J. Khoury, B. Ovrut, P. Steindhardt and N. Turok, Phys. Rev. D 64, 
123522 (2001) 

[5] V. Vanchurin and A. Vilenkin, Phys. Rev.T> 61 (2000) 083507 

[6] N. Chamoun, S. Landau and H. Vucetich, Phys. Lett.B 504, 1 (2001), 
|astro-ph/0008436 

[7] P. de Bernardis et al, Nature 404, 995 (2000). 

[8] S.Hanany et al.,Astrophys. J. 545 (2000) L5. 

[9] H. V. Peiris et al, First Year Wilkinson Microwave Anisotropy Probe 
(W MAP) Observation s: Implications For Inflation, ApJS, 148, (2003) 
213, astro-ph/ 0302225] 

[10] D. Spergel et al, Wilkinson Microwave Anisotropy Probe (WMAP) 
Three Year Results: Implications for Cosmology, |astro-ph/ 0603449 

[11] W. Greiner and A. Schafer, Quantum Chromodynamics, 2nd ed, 
Springer, (1995) 

[12] A.H. Guth, astro-ph/0101507, Proceedings of Cosmic Questions 1999, 
The New York Academy of Sciences Press 

[13] S. Hawking, Phys. Lett. B 115, 4, 295 (1982) 

[14] S. Hawking, Phys. Lett. B 150, 5, 339 (1985) 

[15] R. Brandenberger, ICTP 1997 Summer School in High Energy Physics 
and Cosmology, ed. E. Gava et al, World Scientific, (1997), 
R. Brandenberger, Invited lectures at TASI-94 Colorado, (June 1994), 
|astro-ph/9411049| 

[16] V. Mukhanov, H. Feldman and R. Brandenberger, Phys. Rep. , 215, 
203 (1992) 

[17] Particle Data Group, http://pdg.lbl.gov/ (June 14, 2000) 



11 



[18] C.P.Bachas, Class. Quant. Grav. 17, 951 (2000), |hl^th/0001093[ 

I.Antoniadis an d B.Pioline, Nucl. Phys. B 550, 41 (1999), 
|hep-th/9902055| 



12 



